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Abstract 

We study the electric potential and field produced by disordered distributions of charge to see why 
clumps of charge do not produce large potentials or fields. The question is answered by evaluating 
the probability distribution of the electric potential and field in a totally disordered system that 
is overall electroneutral. An infinite system of point charges is called totally disordered if the 
locations of the points and the values of the charges are random. It is called electroneutral if the 
mean charge is zero. In one dimension, we show that the electric field is always small, of the order 
of the field of a single charge, and the spatial variations in potential are what can be produced by a 
single charge. In two and three dimensions, the electric field in similarly disordered electroneutral 
systems is usually small, with small variations. Interestingly, in two and three dimensional systems, 
the electric potential is usually very large, even though the electric field is not: large amounts of 
energy are needed to put together a typical disordered configuration of charges in two and three 
dimensions, but not in one dimension. If the system is locally electroneutral — as well as globally 
electroneutral — the potential is usually small in all dimensions. The properties considered here arise 
from the superposition of electric fields of quasi-static distributions of charge, as in nonmetallic 
solids or ionic solutions. These properties are found in distributions of charge far from equilibrium. 
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I. INTRODUCTION 



There is no danger of electric shock when handhng a powder of salt or when dipping a 
finger in a salt solution, although these systems have huge numbers of positive and negative 
charges. It seems intuitively obvious that the alternating arrangement of charge in crystalline 
Na^Cl" should produce electric fields that add almost to zero; it also seems obvious that Na^ 
and Cl^ ions will move in solution to minimize their equilibrium free energy and produce 
small electrical potentials. But what about random arrangements of charge that occur in a 
random quasi-static arrangement of charge such as a snapshot of the location of ions in a 
solution? Tiny imbalances in charge distribution produce large potentials, so why doesn't a 
random distribution of charge produce large potentials, particularly if the distribution is not 
at thermodynamic equilibrium? Indeed, some arrangements of charge produce arbitrarily 
large potentials, but as we shall see, these distributions occur rarely enough that the mean 
and variance of stochastic distributions are usually finite and small. More specifically, we 
determine the conditions under which stochastic distributions of fixed charge produce small 
fields. 

The quasi-static arrangements of charge can represent the fixed charge in amorphous 
non-metallic solids or snapshots of charge arrangement of ions in solution, due to their 
random (Brownian) motion. Our analysis does not apply to quantum systems jll], and in 
particular it fails if electrons move in delocalized orbitals, as in metals. Note that the random 
arrangements of charge considered here do not necessarily minimize free energy. 

We consider the field and potential in overall electroneutral random configurations of 
infinitely many point charges. An infinite system of point charges is called totally disordered 
if the locations of the points and the charges are random, and it is called overall electroneutral 
if the mean charge is zero. The configurations of charge may be static or quasi-static, that 
is, time dependent, but varying sufficiently slowly to avoid electromagnetic phenomena: the 
electric potential is described by Coulomb's law alone. In one dimensional systems of this 
type, the potential is usually finite — even though the system usually contains an infinite 
number of positive and negative charges. Even if the system is disordered and spatially 
random, charges of the same sign do not clump together often enough to produce large 
fields or potentials, in one dimensional systems. 

Our approach is stochastic. We ask how disordered can a random electroneutral system 
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be, yet still have a small field or potential. We find the answer by evaluating the probability 
distribution of the electric potential and field of a disordered system of charges. We find that 
the electric field in a totally disordered one dimensional system is small whether the system 
is locally electroneutral or not. The potential behaves differently; it can be arbitrarily large 
in a one dimensional system, but it is usually small in electroneutral systems. 

In two or three dimensional disordered systems, the electric field is not necessarily small. 
We show that in such systems that are also electroneutral the field is usually small. The 
potential, however, is usually large, even if the system is electroneutral. Both potential and 
field are small, if the system is locally — as well as globally — electroneutral (see definition 
below) in one, two and three dimensions. 

We consider several types of random arrays of charges: (a) A lattice with random dis- 
tances between two nearest charges; (b) A lattice (of random or periodic structure) with a 
random distribution of positive and negative charges (charge ±1). Charges in the lattice 
need not alternate between positive and negative, nor need they be periodically distributed; 
(c) A lattice (of random or periodic structure) with random charge strengths. Not all 
charges are ±1, but they are chosen from a set qi, q2, ■ ■ ■ ,qn with probabilities pi,P2, ■ ■ ■ ,Pn, 
respectively, such that 



X^g.Pi = o. (1) 



i=l 



Equation ffl is our definition of electroneutrality in an infinite system. 

n n n 

We use renewal theory [2], perturbation theory p|, and saddle point approximation j^] 
to calculate the electric potential of one dimensional systems of charges and show that it is 
usually small. That is to say, the probability is small that the potential takes on large values. 
Thus, randomly distributed particles produce small potentials even in disordered systems 
in one dimension, if the system is electroneutral. The analysis of one dimensional systems 
requires the calculation of the probability density function (pdf) of weighted independent 
identically distributed (i.i.d.) sums of random variables. This pdf looks like the normal 
distribution near its center, but the tail distribution has the double exponential decay of the 
log-Weibull distribution [5]. We conclude that the electric potential of totally disordered 
electroneutral one dimensional systems is necessarily small, comparable to that of a single 
charge. 

Later in the paper, we define local electroneutrality precisely and show that two and three 



dimensional systems with local electroneutrality usually have small potentials, because the 
potential of a locally neutral system of charges decays like the potential of a point dipole, 
as We show that the potential of typical totally disordered arrays of charges in two 

and three dimensions is infinite even if the system is electroneutral. 

Historically, little attention seems to have been paid to quasi-static random arrangements 
of charge, although much attention has been paid to the equilibrium arrangements 
of mobile charge. In systems of mobile charges, such as liquids and ionic solutions, the de- 
cay of the electric potential may even be exponential, after the mobile charges assume their 
equilibrium distribution. The early theory of Debye-Hiickel 0] shows a nearly exponential 
decay (with distance from a given particle) of the average electric potential at equilibrium, 
originally found by solving the linearized Poisson-Boltzmann equation. In classical physics, 
perfect screening of multipoles (of all orders) occurs in both homogeneous and inhomoge- 
neous systems at equilibrium in the thermodynamic limit, when boundary conditions at 
infinity are chosen to have no effect [l] and there is no flux of any species. This type of 
screeningin electrolytic solutions is produced by the equilibrium configuration of the mobile 
charges 0,0]) which typically takes 100 psec to establish (compared to the 10~^^ time scale 
of most atomic motions) Q| . Many other systems are screened_by mobile charges after they 
assume their equilibrium configuration of lowest free energy such as ionic solutions, 
metals and semiconductors. 

The spatial decay of potential in ionic solutions determines many of the properties of 
ionic solutions and is a striking example of screening or shielding. "Sum rules" of statistical 
mechanics jsi 1^ describe these properties. These rules depend on the system assuming an 
equilibrium distribution, which can only happen if the charges are mobile. 

We consider finite and infinite systems of charges which may or may not be mobile 
and which are not necessarily at equilibrium. We show that the potential of a finite 
disordered locally electroneutral system is attenuated to the potential of a single typical 
charge, whether the potential is evaluated inside or outside a finite system or in an infinite 
system. We note that the behavior of the electric potential and field outside the line or 
plane of the lattice can be analyzed in a straightforward manner by the methods developed 
below. 
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II. A ONE DIMENSIONAL IONIC LATTICE 



Consider a semi-infinite array of alternating electric charges ±g with a distance d between 
neighboring charges. The electric potential $ at a point P, located at a distance R from 
and to the left of the first charge (see Fig. HI) is given by 

~ 4neo\R R + d^ R + 2d R + 3d^' 
g / 1 1 1 



AtteoR V 1 + a l + 2a l + 3a 

AneoR^^l + na ^ ' 

where a = d/R is a dimensionless parameter. The series ^ is conditionally convergent, so 
it can be summed to any value by changing the order of summation jl^. The order of 
summation reflects the order of construction of the system; different orders may 
lead to different potential energies of the system. However, the infinite series 
that determine the electric field are absolutely convergent (see below), so the 
field does not depend on the order of summation of its defining series. Thus, all 
potentials differ from each other by a constant, which presumably refiects the 
different ways the charge distribution could be constructed, while having the 
same electric field. From here on, we consider the ordering in equation (^j) . 

Setting R = d {a = 1) we find the potential at a vacant lattice point (to avoid infinite 
potentials) due to charges located at both directions of the infinite lattice is 

q ^ (-1)""' q 



2$(i? = rf) = 2-^ V = ■ 21og2. 

n=l 

The constant 2 log 2 is known as the Madelung constant of a one dimensional lattice . 

Next we find the asymptotic behavior of the potential $ away from the semi-infinite 
lattice, that is for R ^ d, or equivalently a <^ 1. The following analysis is independent 
of the order of summation of the series (j2I)- Clearly, the infinite sum in eq. ((21) converges, 
because it is an alternating sum with a decaying general term. We expand the potential for 
a -C 1 (away from the lattice) in the asymptotic form 

$ = -^4 (^0 + a\/i + aV2 + ■ ■ ■ ) . (3) 
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The effect of the first charge can be separated from all the others, 

q 1 q I 



where 



Atteq R ineo R + d 
d 



{Vo + dVi + SVa + ■ ■ ■ ) 



(4) 



R+d 1+a 



Comparing eqs.(I2I) and (jH) we obtain 



Vo + aVi + aVs + 



1+a 1+a \l+a 



V2 + ... 



The coefficients Vo,Vi,... are found by equating the coefficients of like powers of a. In 
particular, we find that Vo = 1/2, Vi = 1/4, V2 = 0, so the potential has the asymptotic form 



q 1 

A-keq R 



(5) 



All coefficients Vn can easily be computed in a similar fashion. This result also determines 



the rate at which the potential far away reaches its limiting value, -- 



. The divergent 



2 AnsoR rn 

series for x = 1 has the value Vq = ^ if interpreted as a limit using the Abel sum [l3| 



1-1+1-1+1-1+ 



lim > ( 



lim 



1 



n=0 



»1- 1 + X 



We note that the asymptotic expansion (jS)) can also be found directly from the differential 
equation that the sum 

y{x) = V ^-^x^ 



satisfies 



^-^ 1 + na 

n=0 



axy' + y 



1 



1 + a;' 



(6) 



w.th m.t.al cond.t.on ,(0) ^ 1. The asymptotic form of can easUy be found by standard 
methods In particular. 



lirn y{x) = V - 
x— >i— ' 1 



n=0 



+ na 



The physical interpretation of the asymptotic expansion (0) is that the electric potential 
away from an infinite lattice of charged particles is about the same as if half a single charge 
were located at the origin. The spatial arrangement of the lattice attenuates the effect of its 
charge. The potential near the lattice is determined by a few of the nearest charges and the 



contribution of the remaining charges reduces to that of a half charge placed at a distance 
R ^ d. Obviously, as i? — > the potential becomes infinite, approaching the potential 
produced by just the nearest charge. 

III. ONE DIMENSIONAL RANDOM IONIC LATTICE 

We turn now to solids in which the charges are distributed randomly in several different 
ways. First, consider a semi-infinite lattice of electric charges, in which the sign of each 
charge is determined randomly by a flip of a fair coin. That is, the charges that are located 
at the lattice points Xn {n = 0, 1, 2, . . .) are independent Bernoulli random variables that 
take the values ±1 with probability 1/2. The electric potential of this random lattice is 
given by 

$ 9 ^« 



^ 1 + na ^ 



n=0 

Some discussion of the nature of convergence of the series ((Zj) is needed at this point. The 
convergence of the sum of variances means that the partial sums converge in with respect 
to the probability measure, so the sum ((7j) exists as a random variable $ G L^, whose variance 
is the sum of the variances. Now, the Cauchy-Schwarz i neq uality implies that $ G L^, so 



($) = 0. Note that ((Tj) also converges with probability 1 15 1. 

We use fair coin tossing to maintain the condition of global electroneutrality, though 
arbitrary long runs of positive or negative charges occur in this distribution. Thus some 
realizations of the sequence X„ have runs ('clumps') of substantial net charge and potential. 
The standard deviation of the net charge in a region gives some feel for the size of the 
clumps. The standard deviation in the net charge of a region containing charges is q^/N. 
For large values of A^, substantial regions are not charge neutral. The condition of local 
charge neutrality (defined later) is violated for many of the realizations of charge in this 
distribution. 

Note that a particular set of X„ can produce an infinite potential, despite our general 
conclusions. If, for example, X„ = 1 for all n, the electric potential becomes infinite (see 

oo ^ 

eq. (171)), because > = oo. Nonetheless, the convergence of (171) implies that the 

1 + na 

n=0 

probability that ((Zj) is infinite is 0. In other words, even though the potential is infinite for 
a particular set of X„, the potential is finite with probability 1. This is a striking example 
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of the attenuation of the electric field, even without mobile charge. The attenuation of the 
potential produced by some 'clumpy' configurations of charges occurs even though there is 
no correlation in position, and there is no motion whatsoever. 
The electric field, given by 



E 



n=0 



AttsoR^ ^ (1 + naY 



remains finite for all realizations of because the sum 



s- ^ V L 

AtisoR^ ^ (l + na)2 

n=0 



converges. 

The electric field is bounded (above and below) by 5* and so there is zero probability 
that the function is outside the interval (— S*, S\ The pdf of the electric field is compactly 
supported, even when all charges are positive (or negative). The electric field — unlike the 
potential — is attenuated even if the net charge of the system is not zero, taken as a whole. 

The standard deviation of the field is — ^ — - < \ , which is of the order of 

the field of a single charge at a distance R. 



A. Moments 

The expected value of $ is ($) = 0, as mentioned above. The variance of $ is given by 

-^—\ Y- (8) 

A vacant lattice point in an infinite (not semi-infinite) lattice corresponds to R = d for both 
the charges to the right and to the left. It follows that the variance of the potential there is 
twice that given in (jHJ) with a = 1, that is, 

Varm=2(^\ y\ = 2(^\ ^, (9) 



so that the standard deviation is 
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As expected, the constant / \/?> is larger than the Madelung constant 2 log 2 of the periodic 
lattice, because the potential of the disordered system is larger than that of the ordered one. 

Away from the semi infinite lattice, i.e., for a <^ 1, we can approximate the variance (jH)) 
by the Euler-Maclaurin formula, which replaces the sum by an integral. 



Varm=(^—\ (T- ^-—dx + - + 0(a)] = (^—] (- + - + 0(a] 

so the standard deviation is 



^^\R = ^r^^(^ + 0{a)). (12) 
ATTSoVdR 

The decay law of is more gradual than the decay law — of a single charge. 
Vi? R 



B. The electrical potential as a weighted i.i.d. sum 

The potential (|7j) is a weighted sum of the form (^nXn, where X„ are i.i.d. random 
variables. The distribution of potential is generally not normal. For example, consider 
the weighted sum Yl'^=i 2~"-^n, where Xn are the same Bernoulli random variables. This 
weighted sum represents the uniform distribution in the interval [—1, 1]. It is, in fact equiv- 
alent to the binary representation of real numbers in the interval. Not only does this dis- 
tribution not look like the Gaussian distribution for small deviations, it does not look at all 
Gaussian for large deviations. In fact, this distribution has compact support. It is zero out- 
side a finite interval, without the tails of the better endowed Gaussian. Other unusual limit 
distributions can be easily obtained from sums of the form For example, the weighted 
su m y^, ^^ 3~"X„ is equivalent to the uniform distribution on the Cantor "middle thirds" 
set |l6| in [—1, 1], whose Lebesgue measure (length) is 0. 

Note that the sum 

oo ^ 

(1 + nay+^ 

n=0 ^ ' 

has compact support for every 5 > 0, because the series 

oo ^ 

^(l + na)i+- 

converges for every e > 0. In our case e = 0, so that the limit distribution does not 
necessarily have compact support. Nonetheless, we expect that the probability distribution 
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function of the potential will have tails that decay steeply, even steeper than those of the 
normal distribution. 



C. Large and small potentials. The saddle point approximation 

The existence of the first moment of the sum ((71) depends on its tail distribution, which 

we calculate below by the saddle point method [4]. That is, we calculate the chance of 

finding a pinch of (noncrystalline) salt with a very large potential. For a potential $ defined 

\ -1 



in equation ((7j), we denote the pdf of ( ) $ by f{x). The Fourier transform f[k) 

\A7re0RJ 

of this pdf is given by the infinite product 

fik) = fl cos (j^) , (13) 

which is an entire function in the complex plane, because the general term is 1 + 0(n~^). 
The inverse Fourier transform recovers the pdf 

1 ^ 

fi^) = 7r fik)e^'''dk, (14) 
which we want to evaluate asymptotically for large x. Setting 

g{k, x) = } \ogcos { — )+ikx, (15) 

^ \l + naj 

we write 



00 



/(a;) = ^ y exp{5((fc,x)}rffc. (16) 

The saddle point is the point k for which -^g{k,x) = 0. Differentiating equation p5|) with 

dk 

respect to k, we find that 

k 



00 ^^^^ . 

— c/fc,a; = -> — + tx. (17) 

dk ^-^ 1 + na 

We look for a root of the derivative on the imaginary axis, and substitute k = is. The 
vanishing derivative condition of the saddle point method is then 



tan 



00 tanh 

X = 

n=0 



s 



E Al±^. (18) 



11 



The infinite sum on the right hand side represents a monotone increasing function of s in 
the interval < s < oo, so equation (jlSp has exactly one solution for every x. Near the 
saddle point k = is, we approximate g{k) by its Taylor expansion up to the order 

1 

g{k) ^ g{is) + --^g{is){k - is)\ (19) 
to find the leading order term of the full asymptotic expansion (derivatives of higher order 



of the Taylor expansion can be used to find all terms of the asymptotic expansion ^). We 
use the Cauchy integral formula to calculate our Fourier integral on the line parallel to 
the real k axis through k = is (see Fig. 12)) 

/(a;) ^ i-e^(-) expS^^g"{ts){k-isy^ dk 

= ^e^M / exp <^ g"(ts)— } dz = , (20) 

Equation ()18p has no analytic solution, so we construct asymptotic approximations for large 
and small values of s separately. 



D. Tail asymptotics 

Throughout this subsection we assume that a is small and s is large and we find the tail 
asymptotics of the pdf away from the system (for a ^ 1). For s ^ 1 the Euler-Maclaurin 
sum formula gives 

x= / VI ^»-^/ ^^ + ^tanhg + 0(a). (21) 

g 

Substituting z = , we obtain 

1 + ax 

1 f^tanhz , 1 , , , ^, , 

x = - rf^ + -tanh(s) + 0(a). (22) 

a I Q z 2 



Writing 



tanhz /"^tanhz /""^tanhz — 1 dz 

dz = dz + dz + — 

z Jo z Z Ji z 

, ."^ ianhz , /"°°tanhz — 1 , ^, 

logs+ / dz+ I dz + 0{e~^') 

z .h z 
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we obtain ()22|) in the form 

ax = \ogs + C + ^ + Oia^e-^'), (23) 

where the constant C is given by 

^ tanhz , f°°tanh.z — l , , . 

C = / + / dz. (24) 

Jo ^ Ji ^ 



Exponentiation of equation (|23|) gives the location of the saddle point asymptotically for 
small a and large s as 

g ^ gax-C-a/2+0(a2,e-2»)_ ^25) 

The saddle point approximation (pUj) requires the evaluation of g and its second derivative 
ai k = is. The Euler-Maclaurin sum formula gives 



ai-is) = Vlogcoshj-^ — ) - 

\l + naj 

n=0 ^ ' 



SX 



s log cosh 2; , 1 , , ^/ \ 

- / az -\ — log cosh s — SX + C(as) 

a Jo z^ 2 

s ( /"^ log cosh 2; , r dz /"°° log cosh 2; — 2; , ^/l 
' ' dz^ \ V \ — ^ rf2 + O ' 



Qj V./n .] \ ^ J\ \ S 



H SX + Uias). 

2 2 ^ ^ 

Using equations and (j21|), we find 



g^^s) = C,'^-^-^ + oU^,as], (26) 



where 



^ /"^ log cosh z , /""^ log cosh z — 2; , tanhz , /""^tanhz — 1 , . 
Ci= -^-^ dz+ — dz- dz- dz, (27) 

Jo ^ Ji ^ Jo ^ Jl ^ 

and integration by parts shows that Ci = —1. It follows that 

g{ts) = ---^-^ + o(a,-,as). (28) 

Qj \ d / 
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The second derivative of g is evaluated in a similar fashion 

2 oo 1 - tanh^ ( - — - — 

sr^ \1 +na 



k=is n=0 



'I + nay 



- / (l - tanh^ z) dz (l - tanh^ s) + Oiase'"^') 

as Jq 2 ' 



tanhs . _2, ^_2s^ 



+ 0(ase ,e" 



as 



= -— + 0(as,l,—)e-^'. (29) 
as as 

Substitution of pHjl . (jSHI), and into the saddle point approximation gives 

2vvr 



where the constant C = .8187801402- ■ ■ is given by equation ()24jl . Therefore, the small a 
and large s approximation to the tail of the pdf of $ is given by 

U{x) ~ (31) 



AiieoR / 1 / Aneod /o\ / ^^^o^ /oil 

exp < - I X — C — a/ 2 \ exp < x — C — a/ 2 > > , x ^ oo. 



TT [z\g J a { Q 

It follows from equation (jHlj) that the pdf decays to zero as a double exponential as x — oo, 
which implies that all moments exist. This decay is similar to the extreme value or the log- 
WeibuU (Gumbel) distributions [.5J. The compact support of the distributions of convergent 
series is replaced here with a steep decay. Note also that the decay becomes steeper further 
away from the system, as expected, because the pre-exponential factor of the inner exponent 
is 1/ a = R/d. 

For small x the pdf can be approximated by a zero mean Gaussian with variance Var ($), 
which for small a is 

/.(x)~^^-exp|--^^J |, x-0. (32) 

Near its center, the distribution looks like a Gaussian with a standard deviation that decays 
like 1/\/R, in agreement with equation ()12|). We conclude that the pdf looks normal near 
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its center, but, far away from there, it decays to zero much more steeply, rather hke a cutoff. 
This conclusion is the answer to the question posed in subsection 1111 Bl about the normality 
of weighted sums of i.i.d. random variables. The non Gaussian tails of the distribution are 

IV. RANDOM DISTANCES 

Consider a one- dimensional system of alternating charges without the restriction of equal 
distance between successive charges. In particular, we assume a renewal model, in which 
the distances between two neighboring charges are non-negative i.i.d random variables with 
pdf /(/) and finite expectation value 



^0 

The potential of this random system is also a random variable. 

We show below that away from the system the mean value of the potential V has the 
asymptotic form 



where a = d/R. Equation (jHHj) defines the attenuation produced by the configuration of 
charges. The mean potential of the system is produced by (in effect) half a charge. We 
note that the value 1/2 is exactly the same for both random and non-random systems of 
alternating charges (eq.®). We first note that 



d 





(33) 




(34) 



To find the mean value, we multiply (jH^ by V and integrate (note that < ^ < 
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and then change the order of integration 

(1 



[ f{l) dl h^^oR V Vi{V{R + l) = V} dV 
Jo Jo 



f{l)V{R + l)dl. (35) 
We look for an asymptotic expansion of the form 

V{R) = (Vo + aVi + c?V2 + ...)■ (36) 
Substituting this asymptotic expansion into gives Vq = 1/2 for the 0(1) term, because 

l-a<y^ f{l)^^dl<\. (37) 

The first inequahty is due to the inequahty — > 1 — x. Hence follows. 



V. DIMENSIONS HIGHER THAN ONE 

A. The condition of global electroneutrality 

In dimensions higher than one, global electroneutrality is enough to dramatically atten- 
uate the electric field, but it is not enough to produce a small potential, as shown below. 

Consider the electric potential at a vacant site of random charges located at the points 
of a 2D square lattice 

The variance of $ is 

Var{^) = V , ^ , = 00. (39) 

(n,m)^(0,0) 
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The infinite value of tlie variance means that arbitrarily large potentials can occur with high 
probability. That is, the electric potential is not attenuated. The divergence of the variance 
of the potential of three dimensional systems is even steeper. Therefore, attenuation of the 
potential of totally disordered systems can occur in two or three dimensional systems only 
if some correlation is introduced into the distribution of the locations of the charges. If, for 
example, the signs of all charges alternate, as in a real Na"^Cl~ crystal, the distribution of 
potential will be dramatically different, and greatly attenuated, compared to a two or three 
dimensional system in which many charges of one sign are clumped together. 

The condition of global electroneutrality is enough to ensure the dramatic attenuation of 
the electric field. Indeed, consider a 3D cubic lattice of random charges. The z-component 
of the electric field at a vacant lattice point is 



Xnml COS 

(n,m,07^(0,0,0) 



n 



E = V vVn2 + m2 + /V 



The variance of is finite, 



Var{E,) = Yl 



cos^ 



n 



(n,m,0^(0,0,0) ^ ' 

because convergence is determined by the integral 



< oo. 



2tx I COS 9 sin 9 d9 / — r dr < oo. 



The large potential means that much work has to be done to create the given spatial con- 
figuration of the charges, however, the resulting field remains usually small. 



B. The condition of local electroneutrality 

Here we show that the condition of local electroneutrality implies the attenuation of 
the potential in two and three dimensions. For example, the potential of a two or three 
dimensional lattice of extended dipoles is finite with probability 1, if the orientation of 
dipoles is distributed independently, identically, and uniformly on the unit sphere (see Fig. 

El). 
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Paraphrasing p. 136], we say that a (net) charge distribution p{x) has local charge 
neutrality if the (net) charge inside a sphere of radius R falls with increasing R faster than 
any power, that is, for any x 

lim i?" / p{y) dy = for all n > 0. (41) 

J\x-y\<R 

On a lattice, the number of charges that are assigned to each lattice point can be larger than 
in our example of dipoles (Fig. EJ, thus forming multipoles. The Debye-Hiickel distribution 
also satisfies the local charge neutrality condition. 

The potential of a single lattice point can then be written as an expansion in spherical 
harmonics, if the charges of each multipole are contained in a single lattice box. It can be 
also expanded, if the charge density of each multipole decays sufficiently fast, as 

I 



^OAo(-) = ^ E E ^ ^J^^ (42) 

" /=0 m=-l 

where qim are the multipole moments. In particular, the zeroth order multipole moment is 

q,o = -^J p{y)dy = 0, (43) 



by the condition of local electroneutrality (jHJ: the far potential due to a single lattice 
point decays as (or steeper). The coefficients qim assigned to each lattice point are 
randomized as in the previous sections so their mean value vanishes, meaning that there is 
no preferred orientation in space. (Compare the example of dipoles which do not have a 
preferred orientation.) The mean value of the potential of the entire lattice is then ($) = 0. 
The variance is given by 

yar(<l>) = ^rar($,,,), (44) 

ijk 

where ^ijk is the potential of the charge at lattice point {i,j,k). The potential decays as 
1/r^ (or steeper); therefore the variance decays as 1/r^ = l/(z^ + + /c^)^ (or steeper). The 
convergence of the infinite sum ()44|1 is determined by the convergence of the integral 



1 f°° 1 An 

-,d^ = '^'' -2'^^ = -r<^- (45) 

lr>d ' Jd ' " 

Thus, the variance of the potential is finite and we have shown that local electroneutrality 
produces a dramatic attenuation of potential. As above, the potential away from a charge 
is usually of the order of the potential of a single charge. 
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C. The liquid state 



Screening in the liquid state involves a least three phenomena. (1) The movement of 
charge to a distribution of minimal free energy. (2) The properties of a static charge distri- 
bution with minimal free energy. (3) The properties of any charge distribution. 

If the charge correlation function p{x) minimizes free energy, and is at equilibrium, as 
in ionic solutions, the far field potential is strongly screened. However, the relaxation into 
such a state takes time, typically psec to nsec in an ionic solution under biological conditions 
(see measurements reported in [,10|, and theory summarized in Q]). As long as local charge 
neutrality exists during the relaxation period, the potential changes from attenuated (as 
described above) to exponentially screened, as equilibrium is reached. In fact, the spread of 
potential in ionic solutions has the curious property that it is much less shielded at short 
times than at long times; potentials on the (sub) femtosecond time scale of atomic dynamics 
spread macroscopic distances while potentials on long time scales spread only atomic dis- 
tances. Specifically, potentials on a time scale greater than nano or microseconds spread a 
few Debye lengths, only a nanometer or so under biological conditions, although potentials 
on a femtosecond time scale can spread arbitrarily far depending on the configuration of 
dielectrics at boundaries that govern the violations of local electroneutrality. To make this 
verbal analysis of fast phenomena rigorous, the potentials and fields should be computed 
from Maxwell's equations, not Coulomb's law. 

Non-equilibrium fluctuations may violate local charge neutrality, therefore field fluctua- 
tions can be large. For example, in systems which are not locally electroneutral, potential 
can spread a l ong way, as in the telegraph [2^, Kelvin's transatlantic cable, or the axons 



of nerve cells 



2l|. In such systems, d.c. potential spreads arbitrarily far — kilometers in 



telegraphs; thousands of kilometers in the transatlantic cable; centimeters in a squid nerve 
filled with salt water — even if an abundance of ions (~ 10^^) are present. Local electroneu- 
trality is violated in such systems (at the insulating boundary which separates the inside 
and outside of the cable, e.g., the cell membrane) and that violation allows large far field 
potentials. 
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VI. SUMMARY AND DISCUSSION 



Global electroneutrality ensures the dramatic attenuation of the electric potential and 
field of a one dimensional system of charges. Even if local electroneutrahty is violated, 
and the local net charge is not zero, the potential remains finite in these one dimensional 
systems, even in a random lattice that includes arbitrarily long strings of equal charges. 
Wc have shown that the distribution of the weighted sum of i.i.d. random variables that 
define the one-dimensional electric potential is almost normal near its center, but has very 
steep double exponentially decaying tails. The distances between neighboring charges can 
also be random, without changing the attenuation effect. In higher dimensions, global elec- 
troneutrality is sufficient to dramatically attenuate the electric field, but not the potential. 
However, local electroneutrahty ensures a small potential in two and three dimensions, so 
the electric potential and field is short range in one, two, and three dimensions, if the systems 
are locally electroneutral. 
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FIG. 1: A semi infinite lattice of alternating charges with a distance d between neighboring charges. 
The point P is located at a distance R from and to the left of the first charge. 
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FIG. 2: The integration contour passes through the saddle point k = is in the complex plane. 
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FIG. 3: Two dimensional lattice of dipoles of randomly chosen orientations produce attenuation 
due to the condition of local electroneutrality. 
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